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ABSTRACT 



Context. It has long been suggested that shocks might play an important role in altering the form of the interstellar medium (ISM). 
Shocks enhance gas density and sufficiently dense regions may become self gravitating. Potential star forming clouds within larger 
I molecular clouds, move randomly at supersonic speeds. 

Aims. Depending on the precollision velocity, colliding molecular clouds produce a slab that is either shock compressed or pressure 
confined. In a sequel of two papers (I & II), we simulate molecular cloud collision and investigate the dynamical evolution of such 
slabs. Shocked slabs are susceptible to hydrodynamic instabilities and in the present paper (I) we study the effect of strong shear 
between slab layers on the dynamic evolution of a shock compressed gas slab. Both, head-on and off-centre cloud collisions have 
been examined in this work. We include self gravity in all our simulations. 
Q ■ Methods. Simulations presented here, are performed using the smoothed particle hydrodynamics (SPH) numerical scheme. Individual, 

^ ' pre-collision clouds are modelled as pressure confined Bonnor-Ebert spheres. However, for simplicity the thermodynamic details of 

lyj ' the problem are simplified and the gas temperature is simply evolved by a barytropic equation of state. Obviously, the gas, to some 

, extent suffers from thermal inertial effects. However, we note that the dynamical timescale is much smaller than the local sound 

crossing time so that such effects should have minimum influence. 

Results. Highly supersonic cloud collision produces a cold, roughly isothermal shock compressed gas slab. We find that the shocked 
slab is susceptible to dynamical instabilities like the gravitational instability, Kelvin-Helmholtz (KH) instability and the non-linear 
thin shell instability (NTSI). Growth of instabilities within the slab produces structure in it. The NTSI competes with the gravitational 
instability and the fate of the shocked slab apparently, depends on the relative dominance of either of the two instabilities. Dominance 
of the NTSI causes turbulent mixing between slab layers and dissipates internal energy. Eventually the slab collapses to form a thin 
elongated body, aligned with the collision axis and star formation may commence in it. Our hydrodynamical models discussed here 
suggest that, high velocity cloud collisions may be a viable mechanism for the formation of observed filamentary structure in the ISM. 
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Key words. Molecular clouds (MCs) ~ shock compressed gas slabs - Hydrodynamical instabilities ~ filaments - star formation - 
SPH. 

1 . Introduction Larson (1981) and Elmegreen (2007). Elongated molecular and 

atomic structures are common occurrences in the ISM and a 
Spiral shock waves sweep up matter in galactic arms and pro- number of possible formation mechanisms have been suggested 
. duce large scale structure and dust lanes (Roberts 1969; Fleck (c.f. Kluging & Heiles 1987). One possible mechanism is the 
■ - 1992; Heitsch et al 2006; Elmegreen 2007). Often MCs them- dynamical interaction between fluid flows within MCs (Klessen 
selves are filamentary and such clouds, and may form out of & Burkert 2000; McLow & Klessen 2004). Turbulence decays 
gravitationally unstable gas sheets (Nagai, Inutsuka & Miyama rather quickly, on a timescale equal to the dynamical period of 
1998). Detailed studies of molecular line profiles of MCs have evolution of the cloud. The importance of turbulence in star 
revealed that matter within them is highly unevenly distributed formation has been highlighted by several workers (e.g. Sasao 
and gives them a clumpy appearance, which has led some work- 1973; Klessen, Heitsch & McLow 2000; Klessen & Burkert 
ers to suggest a fractal model of MCs (Elmegreen 1997). 2OOI and references flierein). Gas slabs resulting from coUision 
Schneider & Elmegreen (1977) performed one of the first of stellar winds from binary systems expanding at supersonic 
studies of large scale filamentary structures in the local neigh- speeds become thermally unstable due to rapid cooling, also tur- 
bourhood. Now with improved technology it is possible to map bulent motion within shells makes them unstable to shearing in- 
MCs in greater detail. For instance, dense filamentary structures stabilities (Stevens, Blondin & Pollock 1992). Similarly, power- 
have been observed in the Perseus region where there is ev- ful ionising radiation emitted by young OB associations or su- 
idence for ongoing star formation (HatcheU et al 2005). The pernova remnants sweeps up gas and forms a dense shell that 
Orion Integral Filament (OIF) in the Orion A region has been is gravitationally and thermally unstable (Strickland & Blondin 
studied in great detail and a number of Young Stellar Objects 1995; Dale, Clarke & Bonnell 2007). Low velocity MC coUision 
(YSOs) have been reported to be embedded in it (Chini et al produces gravitationally unstable gas slabs (Bhattal et al 1998). 
1997; Johnstone & Bally 1999). Analogous to random motion of smafler clouds, larger MCs 
Structure within MCs appears to be hierarchical and the (Giant Molecular clouds, GMCs) in galactic arms also exhibit 
densest regions are gravitationally bound, see for instance random motion. Line width of CO emission from GMCs sug- 
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gests that clouds have a random velocity dispersion typically 
of order a few km s"\ superposed on their systematic galactic 
motion (Spitzer 1968; Solomon et al 1987). GMCs are largely 
supported by internal turbulence and have lifetimes of about a 
few tens of Myrs, before being dispersed due to stellar feedback 
(e.g. McLow 1999; Padoan & Nordlund 1999; Klessen, Heitsch 
& McLow 2000; Blitz et al 2007). 

Smaller clouds within GMCs themselves have turbulent ve- 
locities of ~l-2 km s ' on a parsec scale (e.g. Larson 1981; 
Ehnegreen 2000). Collision between such clouds, moving ran- 
domly with supersonic speeds is an efficient mode of dissipating 
kinetic energy. High velocity cloud collision is a violent phe- 
nomenon that might trigger star formation in the resulting com- 
posite object. The gas bodies resulting from dynamically inter- 
acting fluid masses already contain the necessary seeds of dy- 
namic instability. Magnetic field threading the GMCs is gen- 
erally weak, of order a few fiG only, and does not seem to 
play any significant role in damping turbulence within GMCs. 
However, the magnetic field strength could increase dramatically 
due to post-shock enhancement of gas density, which in turn may 
cause further compression of the shocked gas (e.g. Hartmann, 
Ballestros-Paredes & Bergin 2001). 

The problem of colliding clouds has been studied numer- 
ically in the past for instance, by Stone (1970 a,b), Hausman 
(1981), Lattanzio et al (1985) and Hunter et al (1986); more re- 
cently by Bhattal etal (1998) and Klein & Woods (1998), among 
others. Stone (1970 a) studied this problem in one dimension 
by simulating MCs as planes having infinite spatial extent and 
found that post-collision, the clouds coalesced and the result- 
ing gas slab then expanded. The same experiment, in two di- 
mensions also yielded a similar result (Stone 1970b). Hausman 
(1981) performed particle simulations of head-on and off-centre 
cloud colhsions. However, the results were corrupted due to par- 
ticle penetration during colUsion. 

Later Hunter et al, in their grid simulation of colliding flows 
found that such collisions produced slabs which according to 
them were unstable to Rayleigh-Taylor like instability, but could 
not establish the exact nature of this instability. A similar conclu- 
sion was drawn by Stevens, Blondin & Pollock (1995) and Klein 
& Woods (1998). However, both found evidence for non-linear 
growth of instabilities starting in thin shells within the slab. In 
fact the latter suggested occurrence of the non-linear thin shell 
instability (NTSI) in the shocked slab, resulting from cloud col- 
lisions. 

In this paper we demonstrate the occurrence of this insta- 
bility in our SPH simulations of highly supersonic (10 km s ' to 
20 km s"^) head-on cloud collision. For the purpose of this work, 
we use clouds with only token masses, a few orders of magnitude 
smaller than that of real star forming clouds. Although we test 
one model with slightly massive clouds, even their masses are 
still modest. Limited computational resources is a major hurdle 
in using larger number of particles, resolution would be severely 
compromised in case clouds with realistic masses were used. 
The precollision velocity of individual clouds however, is rep- 
resentative of the dispersion velocity of clouds in the galactic 
plane. The possible effect of the SPH artificial viscosity (AV), 
on the growth rate of the NTSI is left for a future paper. 

The plan of the paper is as follows. In §2 we briefly describe 
the SPH code employed for our cloud collision experiments. 
In §3 we discuss the procedure to set up the initial conditions 
for cloud collision experiments and briefly discuss the models 
tested. In §4 we discuss the NTSI and simulation results and 
conclude in §5. Since the column density spans roughly eight 
orders of magnitude, all column density plots presented here are 



plotted on a logarithmic scale. The column density is measured 
in Mopc"^ and the distance in pc. 

2. Smoothed Particle Hydrodynamics (SPH) 

SPH is a Lagrangian interpolation scheme used to treat compli- 
cated problems in astrohydrodynamics. Under the scheme, the 
fluid being investigated is treated as an ensemble of a large num- 
ber of particles or fluid markers. The particles in this case are not 
point masses in the strictest sense, but have a finite spatial extent 
(Gingold & Monaghan 1977; Lucy 1977). 

All the simulations presented here have been performed us- 
ing the DRAGON SPH code. It is a well tested code, the rele- 
vant details of which can be found in (Goodwin, Whitworth & 
Ward-Thompson 2004). It uses the Barnes-Hut tree for gravity 
calculation (0c„y ~ 0.45) and also includes quadrupole moments 
of remote cells. The forces are integrated using a second or- 
der Runge-Kutta integration scheme. Each particle has (50±5%) 
neighbours. The code employs a multiple particle time stepping 
scheme as against global time stepping and uses the standard 
SPH artificial viscosity. 

Extremely dense agglomerations of particles in the compu- 
tational domain are replaced by sinks. A sink is characterised by 
its density, Psm, and radius, Rsm- The concept of a sink particle 
was first introduced in SPH by Bate, Bonnell & Price (1995). 
Rest of the gas particles in the computational domain interact 
with the sink only through gravity. SPH particles coming within 
a predefined radius, called the sink radius (Rsink), are accreted 
by it, provided particles being accreted satisfy certain criteria 
(Bate, Bonnell & Price 1995). In all the simulations presented 
here, the sink density is held fixed at psmk = 10^'^g cm"^. Rsm 
is so chosen that the initial sink mass, Msink, is comparable to the 
minimum resolvable mass in the simulation, M„i„. Now, 

^min ~ ^neibs ' 

where, Nmits is the number of nearest neighbours of an SPH 
particle, m, - ^ and = t^^linkP^'"^- ^^^^ ^® above 
mentioned requirement, 

^ gasPsink ' 

In the present work, the sink particle represents a protostellar 
core. 

3. Numerical experiments: 

3.1. Setting the initial conditions (ICs) 

We first describe the scheme to assemble particles within a MC. 
The MC is modelled as an external pressure confined isothermal 
sphere. An isothermal sphere is characterised by its dimension- 
less radius, ^, and it has been well estabUshed that only those 
spheres having radii smaller than a critical radius, ^cru = 6.45, 
can remain in equilibrium (Chandrasekhar 1939). For our pur- 
poses here, we choose = ^ - which ensures that the sphere 
is supported against self gravity. 

For simpUcity we further assume that the sphere contains 
only molecular hydrogen and if m is the mean molecular mass 
of a gas particle, then the isothermal sound speed ao is defined 
as 

m I 
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Serial 


Experimental 




Pre-collision 


Pext 


Number of 


Head-on 


No. 


details 




Mach number (At) 




particles 




1 


Mmi = = 100 Mo 

Rcldl = Rdd2 = 1 pc, Ti-ij = 


84 K 


35 


Yes 


N,„, = 34893 
yv,,,, = 3516 


Yes 


2 


Mem = Mm2 = 50 Mo 

Rcldl = RcU2 = 0.8 pc, r^-w 


= 54K 


25 


Yes 


N,„ = 129427 
Ngas = 16000 


Yes 


3 


Mcidi = Mcid2 = 400 Mo 




25 


No 


Ngas = 120000 


Yes 




Rcldl = Rdd2 = 2.0 pc, T^id 


= 170 K 








4 


Mcidi = Mm2 = 50 Mo 

Rcldl = Rdd2 = 0.8 pc, Tcld 


= 54K 


25 


No 


f^gas = 40000 


No 
b = 0.2 pc 



Note : The numbers on the left hand side below indicate the respective column numbers of the table. 

(2) Physical details of individual clouds in a simulation. 

(3) The pre-collision Mach number of individual clouds. 

(4) Indicates if clouds are confined by external pressure, Pext, exerted by the ICM particles. 

(5) Total number of particles (Ntot) used in a simulation (ICM -I- GAS). The number of gas particles (Ngas) are separately mentioned. 

(6) Indicates if the clouds were allowed to collide head-on/off-centre; b is the impact parameter of collision. 



Here ks and Tcu respectively, are the Boltzmann constant and 
uniform gas temperature in the sphere. 

The mass interior to this sphere, M{rB), is 



M sphere = Mire) = AnRlpJf^] , 



(2) 



where t//(^) is called the Emden function (Chandrasekhar 1939) 
and the physical radius of the sphere, = Rq^b, Rq having the 
dimension of length is a distance scale factor, and defined as 



Rn 



(4j:Gp,)2 

Taking the ratio of M sphere with we get, 

Msnhere Cl^i^i^B) 



(3) 



' sphere 

rB 



G^B 



2 GM sphered B 

= — — , (4) 

rBHi^B) 
where we define 

The functions if/i^) & //(f) are tabulated by integrating the 
isothermal Lane - Emden equation. The temperature within the 
sphere, Tcu, is calculated using equations (1) and (4). 

Particles within the sphere are randomly assembled using 
three random number generators, each producing a random num- 
ber in the range (0, 1). A random number Hr is chosen such that 



Mir) //(f) 



Msphere I^{^b) 



(5) 



Knowing the RHS of equation (5), the corresponding dimension- 
less radius f can be obtained from the Lane - Emden table. The 
radial distance r of a particle within the sphere is then given by 
r = "Rr^. The probability that a particle will lie in an infinitesimal 
interval of the poloidal angle (0) and {6 + dO) is then, 

sinfl 

pede = —de. 



The total probabihty over the range of the angle 6 is, 



P{0) 



sine' 

fl'=0 ^ 



1 - cos 6 



(6) 



where lig is the second random number Thus, 

cos 6» = 1 - lUe, 

and finally the azimuthal angle, 0, of a particle is fixed as, = 
2TfR^; ^ being a third random number. The (r, 6, 0) coordinates 
are then converted to their cartesian form in the usual way. 
The density of this sphere is 



sphere 



d^ h=^B 



■ e 



(7) 



and ffi = 3, by choice. However, such a sphere in isolation will 
simply diffuse away and therefore needs to be spatially confined 
with a finite external pressure, Pexu acting on its boundary. The 
magnitude of this pressure is 



sphere 



• K^B)e 



-iA(&) 



(8) 



where all symbols have their usual meanings (Chandrasekhar 
1939). 

We simulate the warm external intercloud medium (ICM) 
by assembling particles in an envelope of finite thickness (5/zb). 
Here hs is the average smoothing length of a single SPH particle 
and 5 is a non-zero integer defining the thickness of the envelope. 
Particles in this envelope interact with the gas particles within 
the cloud only by exerting hydrodynamic pressure of magnitude 
Pext, defined by equation (8) above. 

To prevent particles from diffusing across the cloud-ICM in- 
terface, care must be taken to avoid a gradient in particle num- 
ber density. To ensure this, we first calculate the particle number 
density, hjj, near the cloud edge and then evaluate the number of 
particles to be assembled in this envelope, Nenv For the purpose, 
we first obtain the expression for density at the cloud edge in its 
dimensionless form, ps. A trivial mathematical manipulation of 
equation (7) yields. 



Pb 



4;r//(fB) 



(9) 
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Fig. 1. This figure shows the SPH density, ps ph, of the truncated 
isothermal sphere (^b=3). The dashed continuous line is the ex- 
pected density profile, p(?"b), for this sphere, {x = = pspu). 



gas ■ 



1500. 



The particle number density at the edge of the sphere is = 
^ = NgaspB, where m„ the mass of the particle in our test 
sphere is m, - jr— and the volume of the ICM envelope is 
4n[a+s-hB?-i] ^ ^ sphere of unit mass. The number of parti- 
cles to be assembled in the envelope, Ne„v, is then the product of 
this volume with the particle number density, hs- 

The radial coordinate, r, of a particle in the envelope is fixed 
by choosing a random number, "Rr, such that the volume interior 
to this radius is proportional to the volume in the envelope i.e 



'Rr = 



(P - 1) 



iii + d ■ he)^ - ly 



The polar angles i6,(f>) of individual particles are determined 
as before. The polar coordinates of particles are then converted 
to their cartesian form. With this, the assembly of the sphere- 
envelope system is completed. 

To test this scheme of assembUng particles, we assemble 
1500 particles within the sphere (Ngas) and with a choice of 6 
= 5, Nenv - 6159. This assembly is then placed in a periodic box 
and allowed to evolve for a fraction of the sound crossing time. 
The periodic box is simply meant for ghosting particles i.e. par- 
ticles leaving from one face of the cube re-enter from the oppo- 
site face. Fig. 1 shows the SPH density of a truncated isothermal 
sphere, after about half sound crossing time and seems to be in 
reasonable agreement with the expected density profile, marked 
by the dashed continuous curve. 

The box of the ICM particles is set up in a manner similar to 
the envelope of P^xt particles for the truncated isothermal sphere, 
described above. The dimensions of the box are so set that each 
sphere has an envelope of thickness 3hB. Fig. 2 below shows a 
schematic diagram of the initial set up. It shows the precoUision 
MCs barely touching each other in the plane of collision. Circles 
drawn with continuous lines represent the precoUision clouds. 
Each circle is circumscribed by respective dashed circles which 
show the extent of the ICM envelope for each cloud. 

Since we have opted to use physically identical clouds, each 
one of them has the same mass, M^d, radius, Rcm, and an en- 
velope of the ICM particles of thickness dhs. The length of the 
ICM box along the axis of collision (the x axis in Fig. 2) is 

x = 4iR + dhs), 

while that along the y and z axes respectively, is 
y = z^2{R + She), 



and 6 - 3 in the present case. A short coming of the SPH code 
used, is that it requires a cubic periodic box. Each side of the 
cubic box therefore has length x pc. 

Let Ngasi and A^gfls2 be the number of gas particles in each 
cloud (provided externally). The number of particles to be as- 
sembled in the ICM envelope of respective spheres is Nenvi and 
Nenvi- The region marked by stripes in Fig. 2 is that occupied 
by the ICM particles, the volume of which, can be trivially cal- 
culated. We also know the required particle density, n^, in this 
region using which the number of particles to be assembled in 
this region is calculated. The total number of particles in the box 
is simply the sum of Ngasi, Ngasi, Nenvi, Ngmi and the number of 
ICM particles in the striped region of the box. 

Particles in this box are assembled on a HCP lattice. Particles 
within the desired cloud radius are then scooped out and replaced 
with corresponding gas particles. The second cloud is assembled 
in a similar manner. This concludes the process of setting up the 
ICs. 



3.2. Equation of state (EOS) 

In the simulations discussed here, we use a slightly modified 
barytropic EOS given by. 



P(p) 



= (kB/m)x 



( TcA 



gem 



\ lo-wg cm-3 ) 



;p < 10"^^ g cm"^ 

; 10"^^g cm"^ < p < 

5 X 10-2ig cm-3 

;5 X 10"^'g cm-3 <p < 

10-'**g cm-3 

;p > 10"'^g cm"^ 



(10) 



Here, Tcu, is the pre-coUision cloud temperature, Tps, the post 
shock temperature, kg, the Boltzmann constant, m, the mean 
mass of a hydrogen molecule and, y, the adiabatic constant 
which, in the present case is |. This is due to our assumption 
of monatomic gas within individual clouds. 

The density in the thin dissipative region downstream of a 
strong adiabatic shock, is at most four times the preshock den- 
sity, a direct consequence of the Rankine-Hugoniot jump con- 
dition. Shocked gas in astrophysical environs, that can radiate 
efficiently, cools and becomes approximately isothermal. In the 
next few lines we attempt to justify our assumption that the post 
collision gas slab in the models tested here, is approximately 
isothermal. That would be a rather straightforward task in case, 
a post-shock cooling scheme were employed. 

Nonetheless, one may summon for the purpose, the average 
post shock density, temperature and a model for radiative cooling 
proposed by Semadeni et al (2007). The temperature jump con- 
dition, for models 2 and 3 impUes a post shock temperature, Tps, 
of order 6,000 K and 1 1 ,000 K respectively, while that for model 
4 is the same, as in model 2. The cooling rate suggested by the 
latter authors would then correspond to a cooling time, Tcooh 

of 

order a few thousand years only. In that case, the cooling length, 
hooi ikooi ~ Tcooi(io), over which the shocked gas downstream a 
shock cools to the original preshock temperature, would only be 
a tiny fraction of the observed slab thickness in the simulations 
performed. A small cooling length and cooling time (compared 
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Fig. 2. Schematic diagram of the initial set up of colliding clouds. See text above for description. 



to the sound crossing time), presage the assumption of isother- 
mality behind a shock. 

To conclude the argument, let us additionally calculate TcooI 
using the relation, 



2.4 X 10^ 
Tcooi 1 yrs. 



(11) 



Here Hh is the post-shock number density of atomic hydrogen 
(Spitzer 1978). Even by conservative estimates, if we were to 
assume a compression factor of 4, as in a strong adiabatic shock, 
the post shock density would be of order 10^ cm"^, so that equa- 
tion (11) imphes TcooI to be about a few hundred years only, 
much smaller than the local sound crossing time. 

For the models tested here, the post shock density is of or- 
der ~10"^' g cm"-*. Thus, for gas density below this threshold, 
we adopt an isothermal EOS and maintain the gas at the initial 
cloud temperature, Tcu- Following the cloud colhsion, for densi- 
ties above this threshold and less than 5x10"^^ g cm~^, we adopt 
an adiabatic EOS to model the shock. The gas is then held at 
a post shock temperature, Tp,. When the shocked gas becomes 
sufficiently dense (~ Afpdd ~ 10~'^ g cm^^), its temperature is 
brought down to that in the initial, pre-coUision clouds, Tcu. The 
cold, approximately isothermal gas slab is then so maintained up 
to a density of 10"'^ g cm"-'. Above this threshold, the EOS flips 
from being isothermal to adiabatic. This final phase is meant to 
follow putative star forming clumps in the shocked slab. 



3.3. Models tested 

The cloud collision experiments presented here have been listed 
in Table (1) above. The table provides the relevant physical de- 
tails of an experiment. We arbitrarily choose the precollision 
Mach number, M, and the cloud mass, Mcu- The correspond- 
ing cloud radius, Rcu, is determined using the Larson's scaling 
relation between mass and radius. 



Rcidipc) = 0.1 



(PC)( 



Meld 



Mo / 



0.5 



(12) 



(Larson 1981). 

The fate of the post-collision slab depends largely on the ini- 
tial conditions. We observe contrasting features in the evolution 
of the shocked slab in models 1, 2 and 3. While model 1 is un- 
der resolved and just a crude test case, the latter are much better 
resolved, with larger number of gas particles. 

3.3.1. Model 1 

In this case the clouds collide head-on at a highly supersonic 
speed. Post-colhsion the clouds coalesce to form a shock com- 
pressed gas slab, that evolves through the compression and re- 
expansion phases described below. Artificial viscosity dissipates 
the kinetic energy of colliding clouds and heats up the shocked 
slab, whose temperature is denoted as Tps. The gas slab in this 
case does not become Jeans unstable, presumably due to the fol- 
lowing two reasons- (i) insufficient resolution as too few parti- 
cles have been used and (ii) the post-shock density being less 
than that required to support the Jeans instability. 

Consequently, the slab simply re-expands along the collision 
axis. Also, there is no evidence of any other instability in the 
shocked slab. This is obviously due to the poor resolution in this 
model. The fate of such a collision may be envisaged to be a 
diffuse cloud, a part of the ICM. In order to test this outcome, 
we perform a similar experiment in the next case. 

3.3.2. Model 2 

Although the pre-coUision cloud speed in this case is some what 
lower than in the previous case, individual clouds are still highly 
supersonic. Column density plots in Fig. 3 show a time sequence 
of the colliding clouds in this model. Formation of the contact 
layer can be seen in the picture corresponding to t - 0.01 Myr. 
The slab then develops a kink on it (the bending of the slab) 
which then rapidly amplifies, as can be seen in the snapshots 
corresponding to f = 0.03 Myr, t - 0.07 Myr and t = 0.085 Myr. 
Although the slab layers are sheared, the gravitational instability 
dominates in this case. The shocked slab fragments and forms 
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Fig. 3. A column density plot showing a time (in Myr) sequence of the colliding clouds in model 2. The collision produces a shocked 
slab that develops bending modes and kinks in the slab are visible in the plots on the second row. Some material also squirts out, 
jets, from the top and bottom ends of the slab (t = 0.85 Myr). 



a network of filaments. Naively speaking, the slab develops a 
number of holes in it and gas ablates in to filaments. 

Subsequently these filaments become gravitationally unsta- 
ble and produce dense clumps in them. The face-on picture of 
the fragmented slab in Fig. 4 shows filamentary network and 
the clumps within them. The clumps can be identified as bright 
pink blobs in the plot. The gravitational instability grows rapidly 
and the slab fragments on a timescale much smaller than f„ , 
the cloud crushing time defined in §4.2 below. The criterion for 
dominance of the gravitational instability over shearing instabil- 
ities is 



« > 4.7 X 



(13) 



where n is the gas number density in the slab, rj is the ampli- 
tude of perturbation and Antsi is the length of an unstable mode 
(Heitsch, Hartmann & Burkert 2008). 

From the simulation we find that rj ~ 0.02 pc and Amtsi ~ 
0.01 pc, so that equation (13) yields n ~ 10^ cm"^. In the present 
case the post-shock density is ~ 10"'** g cm"^, so that the corre- 
sponding number density is ~ 10^ cm"^, which is nearly an or- 
der of magnitude larger than the calculated n. The shocked slab 
is therefore susceptible to gravitational instability. This outcome 



however, contrasts with that of the previous case where the gas 
slab evolved in to a diffuse cloud. 



3.3.3. Model 3 

This model is similar to the previous two except that the prec- 
ollision clouds here, are unconfined i.e. there are no ICM parti- 
cles. Fig. 5 below shows a time sequence of column density plots 
of the colliding clouds through to the development of the NTSI 
in the shocked slab. The initial behaviour of the post-collision 
gas slab is similar to that in the previous case, eventually the 
NTSI dominates the dynamical evolution of the slab. The grav- 
itational instability thus, seems to be suppressed. Another insta- 
bility to which the slab is susceptible, is the Kelvin-Helmholtz 
(KH) instability. The distinct feature of which is, the formation 
of vortices. Material from colliding clouds streams in from ei- 
ther sides of the contact layer, and smears laterally after striking 
it. Obviously, the gas elements within the slab are in a turbulent 
state and artificial viscosity is responsible for shearing interac- 
tion between them. 

Sure, SPH in its simplest form is unable to intricately re- 
solve the vortex rolls associated with the KH instability but it 
can at least represent the underlying physical process, that of 
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Fig. 4. A column density plot of the fragmented slab in model 
2. We can see a network of filaments with several clumps in 
them, marked as bright pink blobs. The clumps, actually are sink 
particles. Time measured in Myr. 



layers mixing. That, this produces a non-uniform velocity field 
in the plane of the slab is evident from the plot in the top panel 
of Fig. 6. This is a 3-d plot showing the velocity field surface 
in the plane of the slab. Number of secondary peaks on this sur- 
face correspond to the localised circular velocity patterns in the 
slab, obvious features of the KH instability. In the bottom panel, 
the velocity field has been overlayed on the column density plot 
of the slab as seen by an observer stationed in a direction or- 
thogonal to the plane of collision. From this plot we can see that 
the slab develops kinks and the bending modes are more readily 
visible in its central plane. 

With the help of first order perturbative analysis, it can also 
be shown that fluid surfaces in shearing contact develop surface 
gravity waves which give them a rippled appearance also called, 
the breathing modes of the fluid surface. This situation is anal- 
ogous to the formation of ripples on the surface of water bodies 
(Raichaudhari 1999; Shore 2007). Fig. 7 below is a column den- 
sity plot showing a face-on view of the shocked slab. Density 
contours have been overlayed on it to facilitate identification of 
structure within the slab. We can see a dense network of filamen- 
tary structure and clump (p ~ 10"'^ g cm"-') formation is visi- 
ble in regions where filaments intersect. However, this structure 
suffers from tidal disruption due to strong shearing motion be- 
tween slab layers. Klessen & Burkert (2000) and Vietri, Ferrara 
& Miniati (1997) for instance, make similar observations about 
ablation of clumps in their models of turbulent GMC interiors. 

Since we start with unconfined clouds, the shocks confining 
the gas slab move outwards (essentially vacuum in the simula- 
tion) and the slab re-expands and then collapses laterally. See the 
plot of particle positions corresponding to f = 0.50 Myr in Fig. 
12 below. The thermal pressure within the slab is insufficient to 
arrest the collapse, which proceeds towards the formation of a 
filament along the collision axis. Gravitational clustering may 
then commence in this filament. However, we have not followed 
the simulation to that point due to constraints on resources. We 
note that grid codes simulating colliding clouds suggest that the 
re-expansion followed by the collapse phase, is dominated by 
the Rayleigh-Taylor (RT) and KH instabilities (c. f. Miniati et 
al. 1997). However, there is no evidence for the RT instability in 



this model, although vortex formation in the slab suggests exis- 
tence of the latter. 

3.3.4. Model 4 

In this case unconfined clouds collide at a finite impact param- 
eter, b - with other physical parameters the same as those 
in model 2. The supersonic cloud collision results in an oblique 
shock compressed slab. 

Unlike the planar slabs formed in a head-on collision, the 
oblique shocked slab has three distinct regions. The central pla- 
nar region and the two arched extensions at the top and bottom 
ends of the slab. See snapshots in the panels of Fig. 8. These 
extended features are due to the fact that the initial shock front 
effectively blocks the fluid flow from either sides in the central 
regions, while at the top end where the fluid flows from the right 
to left, has higher transverse velocity compared to that in the 
central regions. For the same reason there is also an extension 
at the lower end of the slab, except that the fluid at this end is 
flowing from the left to right. 

Like the former cases, the shocked slab in this case is also 
internally stressed and evolves into a filament aligned with the 
collision axis. Like the previous two cases, the shocked slab in 
this case also quickly develops structure in it. Elongated den- 
sity structures and blobs can be seen in Fig. 9 below, which is a 
column density plot of the slab as seen face-on. 

4. Discussion 

Collision between clouds is characterised by three phases - (a) 
the compression phase, during which post-collision, a contact 
layer forms first and shock waves propagate in individual clouds. 
The slab is confined by ram pressure on its either side which 
strongly compresses the gas within the slab. The thickness of 
the slab depends on the effectiveness of the artificial viscosity. 
Additionally, it could also depend on the efficiency of radiative 
cooling, but such is not the case in the present work as we use a 
plain EOS. 

(b) The next phase is that of re-expansion. During this phase 
the shocks confining the slab propagate outward in the ICM and 
a rarefaction wave propagates back into the slab. This creates a 
central low pressure and low density region in the slab. During 
this phase, some material squirts out from the top and bottom 
of the slab. This is also the phase during which the slab may 
become susceptible to the NTSI and the KH instability due to 
lateral motion of the gas within the slab at supersonic speed (rel- 
ative to the local sound speed in the slab). 

(c) The final phase is that of lateral collapse. The slab expan- 
sion is halted by the ICM by which time, pressure within the slab 
becomes much smaller than the external pressure. Consequently 
the expanding slab material is driven towards the central low 
density region. 

(d) If the externally confining pressure is too weak, the final 
phase could be different and the re-expanding slab may simply 
form a diffuse cloud, amalgamated with the ICM. 

4.1. Instabilities in a shocked slab 

The problem under investigation here is that of the stability 
of a viscous fluid slab confined by accretion shocks on its ei- 
ther sides. Growth of instabilities in a cold, isothermal slab has 
been extensively studied in the past. Linear perturbative analy- 
ses have shown the existence of various unstable modes, see for 
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Fig. 5. A time (measured in Myr) sequence of the column density plots in model 3. Formation of the shocked slab can be seen in 
the snapshot corresponding to f = 0.03 Myr. The NTSI then dominates the slab, as can be seen in the snapshots in the second and 
third rows, after which the slab soon becomes bloated (f = 0. 16 Myr). 

instance Ledoux (1951), Chandrasekhar (1953), Mestesl (1965) 4.1 .1 . The NTSI 



and Simon (1965). However, these investigations assumed slabs 
with no external forces acting on them. 

A shocked slab on the other hand, shows propensity for 
a variety of instabilities. However, the present investigation 
is limited to the gravitational and shearing instabilities only 
since other relevant physical details have been avoided in this 
work. The final state of the shocked slab depends on the relative 
dominance of either of the two instabilities. The shearing 
instabilities are mainly of two varieties- (a) the non-linear 
thin shell instability (NTSI) (b) the Kelvin-Helmholtz (KH) 
instability. The condition for the dominance of the gravitational 
instability over the shearing instabilities, especially the NTSI, is 
given by equation (13) above. 



Growth of instabilities in shocked slabs was first studied by 
Vishniac (1983) by invoking full shock conditions, who sug- 
gested the existence of bending modes of a shocked slab. His 
linear analysis of the growth of these modes showed that the in- 
stability became overstable and led to disintegration of the slab. 
The resulting fragments suffered erosion due to the RT instabil- 
ity. Gravitational instabilities are suppressed in this regime. In 
the limit of a small bending angle, the linear analysis is suffi- 
cient to predict the growth rate of the unstable mode. The lin- 
ear theory however, breaks down when velocities within the slab 
become comparable to the local sound speed and fluctuations in 
column density become of order unity. 

This investigation was extended by Vishniac (1994) to the 
non-linear regime of the bending and breathing modes. The anal- 
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Fig. 6. Top panel : A plot of the velocity field (f = 0.14 Myr) in the plane of the slab in model 3. It is evident that the velocity field 
within the slab is highly uneven and the spikes in it suggest formation of vortices in the slab. Bottom panel : Velocity field in the 
plane of collision overlayed on column density plot of the slab (f = 0.14 Myr) in model 3. The kinks in the slab along with the jets 
are evident from this plot. Velocity is measured in km s"' . 
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Fig. 7. A column density plot of the shocked slab in model 3, as seen in the plane orthogonal to the plane of collision. Density 
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Fig. 8. A time (measured in Myr) sequence of column density plots in model 4. The off-centre cloud collision leads to the formation 
of an oblique shocked slab as can be seen in the pictures of the top row. In the second row we can see that the slab tumbles about 
the z axis and breaks into three parts, with the outer two lobes moving away from the central region. In the bottom row we see that 
the central region evolves into a filament aligned with the collision axis. 
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Fig. 9. A column density plot of the shocked slab in model 4, shortly after its formation. Structure within the slab is evident from 
the plot along with dense clumps, at the intersection of filaments. Time measured in Myr. 



ysis showed that an isothermal shocked slab of inviscid fluid, be- 
came unstable due to non-linear growth of the unstable modes. 
Since this instability commences in thin shells (i.e. between lay- 
ers of the slab), it is called the Thin Shell Instability. The in- 
stability grows rapidly and saturates leading to turbulent mixing 
between layers of the slab. The slab eventually forms a sieve 
like structure and becomes bloated before finally disintegrating 
to form a filament aligned along the collision axis. Gravitational 
clustering can then occur in this filament. A similar observation 
has also been reported by Hueckstaedt (2003) in grid simula- 
tions of colliding fluid flows. According to Burkert & Hartmann 
(2004) a cold, self gravitating gas slab collapses over a timescale, 
tc, and 

^ ^ / Rdd \i 

where Y.cU is the slab surface density. In case of model 3 for 
instance, tc is of order 1 Myr and from Fig. 12 below, we can see 
that the gas slab has collapsed substantially by f = 0.5 Myr. 

Although we start with unperturbed initial conditions, the 
perturbations that develop in the shocked slab are purely from 
white noise. However, the simulation has sufficient resolution to 
resolve the observed NTSI wavelength, A^tsi- In SPH, the spa- 
tial extent of the smallest region resolved is of the order of the 
particle smoothing length, h. An SPH sphere of influence has 
volume lilnh^ 111, so that for a shocked slab with A^^^^ particles, 

where Ax is the shock thickness and remaining symbols have 
their usual meanings. For the chosen dimension of clouds, the 
smoothing length is at least an order of magnitude smaller than 



Ants I- We therefore believe that the features of the instability 
observed in our simulations are physical. 

As a shocked slab accretes material, small perturbations i.e. 
the bending modes of the slab get amplified due to non-linear 
transfer of momentum within the slab. Material streaming in the 
perturbed gas slab deposits a net positive momentum in the con- 
vex regions and a net negative momentum in the concave regions 
of the slab. Transfer of momentum within the slab amplifies its 
bending modes to scales greater than the thickness of the slab 
itself. It is essential to note that the pressure within the slab is 
thermal and isotropic while the external one, is ram pressure. 
This distinction is iiTelevant in the unperturbed state. 

However, once a bending mode of the slab is excited, there 
is a net imbalance in forces acting on the convex and concave 
regions of the slab leading to its amplification. In general, the ac- 
creted material converges in the convex regions of the slab and 
diverges in the concave regions. Slab layers undergo turbulent 
mixing due to strong shear between them. If t is the dynamic 
time scale of the shocked slab (i.e. the time scale on which phys- 
ical properties of the slab evolve dynamically), turbulent mixing 
occurs when t is greater than the time required for generation of 
the next similar layer at a distance greater than the wavelength 
and having wavenumber k. If L is the thickness of the slab and t 
is its age then 

T » (14) 
kL 

is the condition under which turbulent mixing is suppressed 
Vishniac (1994). 

If T ~ f, the turbulent mixing is weak and is only marginally 
effective in the post-shock slab. The turbulent mixing of slab lay- 
ers gives rise to eddies on the scale of slab thickness which is an 
effective means of energy dissipation. The NTSI grows at a rate, 

~ aokikxj'^ ■ C-/'^ (15) 
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Vishniac (1994). 

Here Xm denotes the mid-plane of the slab and Ca is a con- 
stant of order unity. This constant describes the extent of dissipa- 
tion resulting from mixing between slab layers. Klein & Woods 
(1998) estimate 2 < Q < 4 and a rate of saturation > aokikL)^^^. 

4.2. Observations from simulations 

Clouds colUde and first form a layer of contact after which, 
shock waves propagate into individual clouds. The cloud ma- 
terial streams into the layer and the timescale on which a cloud 
is destroyed is called the cloud crushing time, r^, . Shock waves 
propagating in a cloud reach its far edge in this much time, and 
tcr ~ ^ (Klein & Woods 1998). Here v, is the velocity of the 
shock wave. 

The density and pressure (thermal + ram) gradients respec- 
tively, have opposite signs across the shocked slab. The slab sur- 
faces are shocked while pressure within the slab is thermal. Thus 
the pressure gradient across the slab surface has opposite sign 
relative to the density gradient, that increases steeply. The shocks 
confining the slab apparently arrest the inflow of gas from the in- 
dividual clouds. The colliding clouds are therefore held up. 

This can be seen for instance, in model 2 (§3.3.2). The col- 
umn density plot corresponding to t = 0.085 Myr (which is 
greater than tcr) in Fig. 3 above, shows that the colliding clouds 
survive even beyond their crushing time, tcr- In fact at this epoch, 
the shocked slab has undergone gravitational fragmentation. Fig. 
4 shows the column density plot (f = 0.085 Myr) of the frag- 
mented slab. Filamentary structure within the slab is evident in 
this plot. The filaments then develop condensations and form a 
number of clumps in them. These clumps can be seen as bright 
pink spots in the picture. 

The clumps are actually sink particles and we remind our 
reader that sinks, in the present work represent protostellar cores. 
Soon after its commencement, core formation spreads like an 
epidemic and within a few thousand years, sixteen cores are 
formed in various filaments, whence the simulation was termi- 
nated. Cores accrete material channelled along their natal fila- 
ments. 

Immediately after its formation, the shocked slab develops 
kinks and features of the NTSl are visible in it. The bending 
modes of the slab subsequently amplify and soon, it becomes 
bloated. This can seen in the snapshots corresponding to t = 0.03 
Myr, t = 0.05 Myr and t = 0.07 Myr of Fig. 3. Interestingly, there 
is a confluence of the gravitational instability and the NTSl but 
the former dominates over the latter and the condition for dom- 
inance of the gravitational instabihty is given by equation (13) 
above. Thus, we note that the gravitational instabihty can domi- 
nate the shocked slab even in presence of the NTSl, provided the 
slab is sufficiently dense. 

Although we have seen evidence for clump formation in this 
model, only the self gravitating clumps survive. Other smaller 
clumps either merge to form larger clumps or simply diffuse 
away. We therefore envisage this model as one leading to the 
formation of filamentary structure, that may further spawn star 
formation. This is therefore an interesting case from the perspec- 
tive of star formation. However, we do not observe any evidence 
for secondary fragmentation of cores in this model. 

In model 3 we perform a similar experiment with higher 
resolution but no intercloud medium. The shocks confining the 
post-collision slab move outward (into vacuum) and thus, permit 
re-expansion of the slab, as can be seen in the pressure-distance 
plots of Fig. 11. Distance in parsecs is marked on the x - axis 



and pressure, normalised relative to the initial pressure, has been 
marked along the y - axis. The plots in the top row (/ = 0.48 
tcr), correspond to the post collision jump in density and pres- 
sure, respectively. The plot of pressure shows the formation of 
shocks at the slab surfaces while that within the slab is thermal 
and uniform. 

The pressure and density gradients respectively, have oppo- 
site signs across the slab surfaces, making it unstable. Unlike 
Hunter et al (1986), we are however, disinclined to term this sit- 
uation as Rayleigh-Taylor (RT) unstable. The outwardly propa- 
gated and weakened shocks can be seen in the plots correspond- 
ing to ? = 1 1 .4 ?cr- During this phase the central density and pres- 
sure within the slab falls, as can be seen in the plots in the second 
and third row of this figure. Spikes in the density distribution 
indicate clump formation which obviously, do not survive. The 
slab continues to re-expand and flattens laterally. Eventually it 
becomes elongated with a large aspect ratio and its long axis 
lying in the plane of collision. Fig. 12 shows a plot of particle 
positions in the plane of collision. 

As the shocked slab evolves, a few breakaway blobs are vis- 
ible near its surface. These can be seen in the column density 
snapshot corresponding io t - 0.16 Myr in Fig. 5. This is per- 
haps due to the inability of SPH in its simplest form to model the 
KH instability, apparently due to formation of artificial bound- 
aries between regions having steep density contrast (Agertz et 
al 2007). This problem can be fixed by including some addi- 
tional correction terms to account for thermal energy dissipation 
across these boundaries (Price 2007 11). Our code does not in- 
clude these corrections, however this shortcoming does not seem 
to have any adverse bearing on the results of the simulations and 
in any case it is not our aim here, to model the KH instabilities 
in their minute details. 

Similar breakaway blobs, having much smaller densities (at 
least by two orders of magnitude), as compared to those ob- 
served in our simulations, have been reported by Heitsch et al 
(2008). They interpret this observation as evidence for core for- 
mation. While these latter authors began with perturbations on 
lateral faces of colliding fluid streams and grid codes apparently, 
being able to better handle density contrasts, it is plausible that 
their observation of core formation is physical however, we pre- 
fer to discount the blobs as numerical artifacts despite they hav- 
ing densities about an order of magnitude larger than that of the 
slab. 

Further, as in the previous case, the shocked slab in this case 
also shows features of the NTSL See column density plots corre- 
sponding to f = 0.1 1 through to 0.16 Myr, in Fig. 5 above. These 
pictures show that the kinks (bending modes) in the slab, grow 
rapidly and soon it becomes bloated (f = 0.16 Myr). In the bot- 
tom panel of Fig. 6, the local velocity field has been overlayed 
on the column density plot of the shocked slab (f = 0.14 Myr), 
by which time features of the NTSl have developed well. This 
picture shows amplified kinks in the slab. By measuring the am- 
plitude of the perturbation and the slab thickness, we calculated 
the growth rate of the NTSl as 

~ (0.56)(feL)°•■'^ (16) 

where k is the wavenumber of the unstable mode and k = -r^, 
L is the thickness of the slab. Fig. 10 shows the growth rate of the 
NTSl. The measured growth rate has been plotted in red while 
the smooth green curve is a power-law fit (equation 16) to the 
one in red. 

The observed growth rate of the NTSl in this simulation is in 

good agreement with that reported by Blondin & Marks (1996) 
and Klein & Woods (1998). They suggest a growth rate of order 
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Fig. 10. The growth rate of the NTSI in model 3. The red curve 
corresponds to the actual growth rate of NTSI in the simulation 
while the green curve is a power law fit (Eqn. (16)) to the ob- 
served growth rate. Time on the x axis is measured in units of 
the cloud crushing time t^r- 

{kLf ''. These respective growth rates of the NTSI however, do 
not match with that proposed by Vishniac (1994) and given by 
equation (15) above. This later equation predicts a much faster 
growth (A:'^), of the instability. In fact the instability in this case, 
has been observed to grow on a timescale comparable to the 
sound crossing time across the shocked slab. 

We have argued that gas within the shocked slab is in turbu- 
lent motion, so that the slab is internally shocked and mechanical 
energy dissipated. Accompanied with this, is the loss of linear 
momentum of gas elements within the slab. The larger the vis- 
cosity, the greater is the dissipation. We have already noted that 
growth of the NTSI is related to momentum transfer between 
perturbed regions of the slab. This dissipation should heat the 
gas up, a phenomenon that is obviated due to a bary tropic EOS. 
This is followed by lateral expansion of the gas slab, which soon 
becomes bloated. Rounding of the perturbed surfaces will reduce 
local pressure imbalances, another crucial factor inconducive to 
the growth of NTSI. Concurrence of these two factors may be 
responsible for arresting the growth of this instability. The NTSI 
is observed in models 2 and 3 and grows at a comparable rate in 
either cases. Artificial viscosity being the only common factor 
between the two cases, is therefore likely to have influenced its 
growth. Other cited literature here however, does not discuss the 
possible reasons for the observed slow growth rate. Looking at 
the end product of this simulation, it may be suggested that such 
a collision is likely to produce a diffuse cloud which may col- 
lapse further and form a long filament. Gravitational clustering 
may commence in this filament. 

4.2.1. Off-centre cloud collision 

The shocked slab is internally stressed and thus exhibits a non- 
uniform velocity field. Internal shear renders the slab layers KH 
unstable. The picture in Fig. 9 is a column density plot of the 
slab soon after its formation. The slab forms well defined struc- 
ture which is either filamentary or globular. Some small globules 
merge with other similar globules to form larger ones while a 
few others are destroyed. As in model 3, the colliding clouds and 
therefore the post-collision gas slab, are unconfined. The confin- 
ing shocks can therefore propagate outward. Thus, the slab in 
this case as well, passes through the same phases as the one in 
the previous case. Additionally, it also exhibits rotational mo- 
tion about an axis orthogonal to the plane of collision iz axis). 



As the slab tumbles, two satellite lobes break off and move away 
radially. See the snapshot corresponding to f = 0.27 Myr in Fig. 
8. 

The slab dissipates thermal support through internal shearing 
motion. Usami, Hanawa & Fujimoto (1995) provide a detailed 
analytic account of oblique shock compressed slabs, according 
to which gravitational instabilities are suppressed in such slabs. 
Shearing motion increases the effective sound speed and there- 
fore damps the growth rate of gravitational instability, which is 
probably exacerbated by the type of EOS employed in this work. 

Further, these authors also predict the end product of such 
slabs to be an elongated object, similar to the one we observe 
at the end of this calculation. We terminated this calculation af- 
ter the formation of a filament in the plane of collision. The se- 
quence of colliding clouds and then formation of the filament has 
been shown in a time sequence of column density plots in Fig. 
8 above. We note that the slab in this case evolves on a much 
longer timescale as compared to the earlier models. This is pre- 
sumably due to the additional support provided by the angular 
momentum associated with the slab rotation. 

We notice some perturbations on the slab surface (see snap- 
shots for f - 0.07 Myr and f = 0.12 Myr in Fig. 8) but do not 
amplify as much as those in the earlier model. The NTSI thus, 
seems to be damped out. The Coriolis force acting on the rotat- 
ing slab may suppress perturbations normal to the slab surface, 
in other words, the bending modes of the slab. This is however, 
contrary to the findings of Whitworth et al (1995) who simu- 
lated off-centre low velocity cloud collision and found that the 
pressure compressed oblique slab, underwent gravitational frag- 
mentation. However, with just a few thousand particles it is un- 
likely that they would have been able to resolve the gravitational 
instability. 



5. Conclusions 

Cloud collision plays an important role in the interstellar gas 
dynamics and is one of the mechanisms of dissipating energy. 
Sufficiently dense clumps resulting from fragmentation of the 
post collision gas slabs become self gravitating and spawn sin- 
gle/multiple star formation. A violent phenomenon like cloud 
collision also replenishes the pool of elements in the galac- 
tic disk. Post shock gas attains incredibly high temperatures 
(> lO"* K) however, if the shocked gas radiates efficiently, it 
cools rapidly and the temperature soon drops, sometimes even 
to the initial preshock temperature. Gas is further compressed 
by the shock and postshock gas bodies may attain densities as 
high as ~ 10'*cm"'' - 10^ cm""*. Under favourable conditions, 
star formation may commence in these gas bodies . Thus, gas in 
the ISM goes through a cycle of phase transition starting from- 
diffuse gas to MCs to shocked slabs to star forming cores to stars 
and back to diffuse gas. In some cases however, stars just may 
not form between the shock phase and the final diffused state. 

The scope of the cloud collision problem is quite large and 
here we attempt to investigate only its dynamical features. With 
the aid of the simple models tested in this work, we have at- 
tempted to suggest that cloud collision provides a viable mech- 
anism for the formation of elongated structure in the ISM. The 
densest regions in these filaments may become self gravitating. 
We have concentrated only on the gravothermal aspect of the 
matter, thereby discounting the effect of the interstellar chem- 
istry. We have tried to ascertain the role of dynamical instabili- 
ties on the evolution of shocked gas slabs and have shown that, 
there is much more to the problem than just a balance between 
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Fig. 11. Time evolution of density and pressure against distance along the collision axis in model 3. Time is measured in units of 
the cloud crushing time, fee Both, density and pressure have been normalised relative to their respective initial values. In the top 
panel we can see the jump in density and pressure, immediately after cloud collision. Note that the average density contrast behind 
the shock is roughly four, as should be the case for a strong adiabatic shock. The density behind the shock increases further, as we 
move in to the isothermal shock approximation, see Eqn. 10 above. The shocks then move outward and the slab begins to re-expand, 
accompanied with a drop in density and pressure within the slab. We can see this in respectively, the second and third panels above. 



gravity and thermal pressure. Fig. 13 below is a plot of the nor- 
malised gravitational and thermal energy against time, for mod- 
els 2, 3 and 4. A common feature of these plots is that, post- 
shock, the thermal energy rises steeply, attains a maximum be- 
fore rapidly falling off during the re-expansion phase of the slab. 
However, the shocked slab in model 2 expands only marginally, 
and undergoes gravitational fragmentation. The slab in this case 
is thermally supported, as is evident from the red curve in Fig. 
13. The respective slabs in models 3 and 4, re-expand consider- 
ably and eventually end up at a temperature close to that in the 
precoUision clouds. The gravitational energy, in models 3 and 4, 
is almost in sync with the thermal energy; with model 2 being 
an exception. The cold shocked slab in this latter case becomes 
gravitationally unstable and fragments. The steep rise in gravita- 
tional energy corresponds to self gravitating fragments. 



In all the models tested here, the shocked slab becomes 
bloated soon after its formation which is a unique feature of 
the NTSl. The NTSI along with the KH instability, play a vital 
role in energy dissipation via vortex formation within slab lay- 
ers. Evidence for the NTSl has also been observed in grid sim- 
ulations of colliding turbulent fluid flows (Heitsch et al 2008). 
However, the gravitational instability dominates a shocked slab 
provided it is sufficiently dense else it is suppressed and the 
slab re-expands. This can be seen in models 1, 3 and 4. Usami, 
Hanawa & Fujimoto (1995) for instance, make a similar claim. 
The NTSl appears to dominate the unconfined slab in model 3, 
while model 1 is poorly resolved, although the final result is the 
same. The slab in model 2, contained by the ICM, is more inter- 
esting as we observe the NTSl competing with the gravitational 
instability and finally the slab fragments. The oblique shocked 
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Fig. 12. A plot of particle positions showing the formation of an 
elongated object in model 3. At this epoch (f = 0.5 Myr), the 
shocked slab has re-expanded and undergone significant flatten- 
ing. 

slab in model 4 develops kinks which seem to be damped and 
so, nothing conclusive can be said about growth of the NTSI in 
this model. 

Klein & Woods (1998) and Heitsch et al (2006) also report 
similar conclusions. The former additionally note that, the NTSI 
is suppressed in adiabatic shocked slabs as such shocks do not 
radiate post-shock kinetic energy and therefore tend to rapidly 
re-expand, on a timescale much shorter than that required for 
the growth of bending modes. We measure the growth rate of 
the NTSI in one of our models, which agrees very well with that 
reported in literature elsewhere. 

In the present work, we have adopted a simple EOS to ac- 
count for radiative cooling and the assumption of an isother- 
mal shock may have led to excessive post-shock compression. 
It is therefore essential to employ better cooling techniques (e.g. 
Koyama & Inutsuka 2000; Vazquez-Semadini et al 2007). It is 
also essential to investigate the effect of SPH artificial viscos- 
ity on the growth rate of NTSI. We adopted the standard choice 
of viscous parameters in all our simulations here. The models 
presented here are extremely simple, since we have preferred 
to use physically identical clouds without any internal perturba- 
tions. Also, we have not included magnetic fields in the present 
work. It would be an instructive exercise to start with magne- 
tised clouds and then collide them. We can then study the effect 
of magnetic field on the evolution of the shocked slab and in 
particular, on the growth of various dynamical instabihties in a 
magnetised gas slab. 
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